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Dilepton rapidity distribution in the Drell-Yan process at NNLO in QCD 



m 
o 
o 

(N 

C 
3 



> : 

<N . 
G\ ■ 

\o : 
o . 

m ■ 

o ; 

^ : 

a: 
i 

^ ■ 



X 



Charalampos AnastasiouJ Lance Dixon, 1 Kirill Melnikov, 2 and Frank Petriello 1 

1 Stanford Linear Accelerator Center, Stanford University, Stanford, CA 94309, U.S.A. 
2 Department of Physics & Astronomy University of Hawaii, Honolulu, HI 96822, U.S.A. 

We compute the rapidity distribution of the virtual photon produced in the Drell-Yan process 
through next-to-next-to-leading order in perturbative QCD. We introduce a powerful new method 
for calculating differential distributions in hard scattering processes. This method is based upon 
a generalization of the optical theorem; it allows the integration-by-parts technology developed for 
multi-loop diagrams to be applied to non-inclusive phase-space integrals, and permits a high degree 
of automation. We apply our results to the analysis of fixed target experiments. 
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The production of lepton pairs in hadronic collisions, 
known as the Drell-Yan (DY) process Q, was the first 
application of parton model ideas beyond deep inelas- 
tic scattering. Due to its clean theoretical interpretation 
in terms of quark-antiquark annihilation into a vector 
boson, and large event rates, the DY process has been 
studied extensively, and will continue to be investigated 
at both the Tevatron and the LHC. The DY process 
provides valuable information about parton distribution 
functions (pdfs), enables measurements of the produc- 
tion rates and masses of W and Z bosons, and furnishes 
a sensitive test for many varieties of new physics, such as 
the additional gauge bosons that appear in many exten- 
sions of the Standard Model. It will also be used for the 
more prosaic purpose of monitoring partonic luminosities 
at the LHC. 

Despite the importance of the DY process and the 
significant amount of work devoted to its description, 
the calculation of higher order QCD corrections has pro- 
ceeded slowly. The next-to-leading order (NLO) QCD 
corrections to the total cross section, and the xf and 
rapidity distributions, were calculated nearly 25 years 
ago 0; the NNLO corrections to the total cross section 
were obtained eleven years later Q. No complete calcu- 
lation of the NNLO QCD corrections to any differential 
distribution has been performed, although partial results 
exist 0. 

Recently, the NNLO virtual corrections to several in- 
teresting hard scattering processes in QCD have been 
computed |5j; however, the calculation of real emission 
contributions, required for complete NNLO predictions, 
is still in progress. These contributions entail a care- 
ful analysis of perturbative multiparticle final states in 
generic hard scattering events. While it is certainly use- 
ful to solve this problem in complete generality, it is also 
useful to study specific examples, especially those most 
urgently needed in experimental analyses. It is possi- 
ble to develop alternative methods of calculation which 
can be used to compute basic differential distributions. 



In Refs. |E it was shown how to combine the opti- 
cal theorem with multi-loop computational methods to 
compute phase space integrals. In this Letter we present 
a non-trivial application of these ideas; we compute the 
rapidity distribution of the virtual photon produced in 
the DY process through NNLO in perturbative QCD. 

We shall apply our results to the production of lep- 
ton pairs in proton-proton (pp) collisions at center-of- 
mass energies accessible to fixed target experiments. The 
most recent measurements come from Fermilab experi- 
ment E866/NuSea, which measured the dimuon produc- 
tion cross section in pp and proton-deuteron collisions at 
i/s « 40 GeV for muon invariant masses in the inter- 
val 4—16 GeV [g. These experiments are sensitive to 
both the x — ► 1 components of the valence quark distri- 
bution functions and the moderate x components of the 
sea quark distribution functions of the proton. Neither 
of these kinematic configurations are well constrained by 
other data, so the E866 measurements provide valuable 
input to a global pdf fit. The precision of the E866 mea- 
surement is better than 10% per bin. Given the sig- 
nificant (~ 40%) NLO corrections at such energies, the 
complete NNLO computation is required. Although in 
principle both photon and Z boson exchanges contribute 
to this process, the Z exchange component is suppressed 
by M 2 /M§, where M is the invariant mass of the lepton 
pair. This effect is approximately 1% for the relevant 
invariant masses, and will be neglected in our analysis. 

The NNLO calculation is quite challenging technically. 
Existing techniques for computing phase-space integrals 
are incapable of handling problems of this complexity. 
We introduce here a powerful new method: we extend the 
optical theorem in such a way that the calculation of dif- 
ferential distributions becomes possible using techniques 
developed for multi-loop calculations. To achieve this, we 
represent the rapidity constraint by an effective "prop- 
agator." This propagator is constructed so that when 
the imaginary part of the forward scattering amplitude 
is computed using the optical theorem, the "mass-shell" 
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constraint for the "particle" described by this propaga- 
tor is equivalent to the rapidity constraint in the phase 
space integration. We then use the methods described 
in Ref. for computing inclusive cross sections, keeping 
the fake particle propagator in the loop integrals, and 
deriving the rapidity distribution as the imaginary part 
of the forward scattering amplitude. 

The production of a lepton pair in a high-energy 
hadronic collision occurs in two distinct steps: first the 
quarks and gluons from the colliding hadrons annihi- 
late to create a highly virtual time-like photon; then the 
photon decays into a pair of leptons. In the center-of- 
mass frame, the two colliding hadrons have momenta 
-Pi, 2 = \/s/2 (l,Oj_,±l). A virtual photon of invari- 
ant mass M produced in the collision has momentum 
P-y — (E,p±,p z ). Its energy and momentum are related 
by the "mass-shell" condition E 2 — p 2 — p 2 = M 2 , while 

its rapidity is defined as Y = — In 
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We first compute the partonic hard scattering cross 
sections, and then convolute them with the pdfs of the 
colliding hadrons. The partonic rapidity distributions 
for the hard scattering of partons with momentum 
Pi — x\P\ and p2 — X2P2 respectively, are obtained by 
integrating the hard scattering matrix elements over the 
phase-space of the final state particles with the rapidity 
and mass of the virtual photon kept fixed: 
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The rapidity constraint can be rewritten using the in- 
coming parton momenta as 
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The simple Lorentz boost properties 



where u = —e 2Y . 
of the rapidity are helpful at NNLO. In comparison to 
the total cross section computation, only one additional 
dimensionless variable, it, is introduced. Computation 
of the distribution in xf = 2p z /y/s, for example, would 
require two new variables. 

At leading order in a s , the production of the virtual 
photon occurs through the annihilation of a qq pair. Only 
the virtual photon is produced in the collision, rendering 
the phase-space integrations trivial. At higher orders in 
a s , inelastic channels contribute. At C(a s ), for example, 
we must consider also qq — > 7*5 and qg — > 57*. It is still 
quite simple to perform these phase space integrations 
using standard techniques. At higher orders, however, 
this approach becomes impractical, so we adopt instead 
the method of Ref. 0, which can be applied efficiently 
at NNLO. 

We first represent the <5-function in Eq. @ as the imag- 
inary part of an effective propagator: 
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Next we map the constrained phase-space integrals onto 
forward scattering loop integrals |g. We denote the dif- 
ference of propagators with opposite iO prescription, such 
as that shown in Eq. @, by a cut propagator; final state 
particles on mass-shell are also represented by cut prop- 
agators. The 5- function constraint becomes an un- 
conventional propagator, linear in the loop momentum. 

At NLO, we must consider integrals of the following 
general form: 



I(Vl,U 2 , ...,V 5 ) = 
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where A t = k 2 - M 2 ± iO, A 2 = (k+p 1 ) 2 , A 3 = (k + 
Pi+P2) 2 ±«0, At = (k+p 2 ) 2 , and A 5 = k-pi~uk-p 2 ±iO. 
The propagators Ai,A 3 and A 5 should be "cut" accord- 
ing to Eq. (PJ), indicating that the corresponding par- 
ticles are on-shell. The propagators A1..5 are linearly 
dependent; we can therefore eliminate both A2 and A4 
from the integrand in Eq. by partial fractioning. Par- 
tial fractioning produces integrals with either v\ , 1/3 or 1/5 
equal to zero. When the cutting rule © is applied, such 
integrals vanish. All phase-space integrals of the form 
of Eq. @ can be reduced to a single "master" integral, 
7(1,0,1,0,1). The fact that only partial fractioning rela- 
tions are required to perform this reduction is specific to 
NLO; we will discuss a more general reduction technique 
when we consider the NNLO corrections. 

We compute the virtual corrections to the leading or- 
der production process qq — ► 7* in the standard fashion, 
since the rapidity constraint leaves this calculation unaf- 
fected. After combining the real and virtual corrections 
and performing the collinear factorization, we arrive at 
the LO and NLO results for the partonic rapidity distri- 
butions 0], which we present here for completeness. 

We write the partonic differential cross section for the 
process i + j — > 7* A, renormalized in the MS scheme, as 
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where a s = a s {M) is the strong coupling constant as- 
suming rif massless quark flavors, renormalized at the 
scale M. The factorization scale is also set equal to M; 
the dependence on both the renormalization and the fac- 
torization scales can be restored by using the renormal- 
ization group invariance of the hadronic cross section. 

At the lowest order in a s , the virtual photon can be 
produced only in the collision of a quark and antiquark 
of the same flavor. Therefore, 
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where z = M 2 / s, s is the partonic Mandelstam invariant 
and y = (u — z)/(l — z)/(l + u). 

At NLO, the qq channel receives 0(a s ) corrections, 
and the qg and qg channels contribute. We find 
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Results for the other channels arc obtained by permuting 
partonic labels and changing y — ► 1 — y in Eqs. 

We now discuss the calculation of the NNLO contribu- 
tions. The purely virtual correction to the rapidity distri- 
bution is the same as the virtual correction to the total 
cross section, and is straightforward to compute using 
standard techniques. We compute both the real- virtual 
and the real-real corrections using the method described 
above. However, to achieve a complete reduction to 
master integrals at NNLO we must supplement the par- 
tial fractioning identities with integration-by-parts (IBP) 
identities typically used in the reduction of loop in- 
tegrals. Our substitution of the rapidity constraint with 
an effective propagator facilitates the use of IBP iden- 
tities in phase-space integrals. The combined system of 
equations reduces all required integrals to approximately 
thirty master integrals, which depend on two variables, u 
and z. The equations are also used to construct differen- 
tial equations satisfied by the master integrals, following 
Refs. |(| [^}. Here, two first order inhomogeneous par- 
tial differential equations are derived for each integral. 
These equations are solved, and the boundary conditions 
obtained by considering simple kinematic limits. 

At NNLO, the following partonic channels contribute: 
qq, the scattering of a quark and anti-quark of the same 
flavor; q(q)g; gg; and ?i (fa (<£,), the scattering of quarks 
(anti-quarks) of arbitrary flavor. The complete analytic 
results for the partonic cross sections are quite lengthy, 
and will be presented elsewhere. 

Integrating the partonic cross sections over the virtual 
photon rapidity, we reproduce the 0(0?) correction to 
the total cross section computed in Ref. {3| . This provides 
a strong check on our result. 

Finally, we must convolute the renormalized partonic 
cross sections with the partonic structure functions to 
obtain the experimentally measured cross section. We 
consider the doubly differential cross section d 2 cr/dAfdY": 
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where z = M 2 /{sxxx 2 ), u — (xi/x 2 )e~ 2Y , and a is the 
electromagnetic coupling evaluated at the scale M; nu- 
merically, or 1 w 132. We use a consistent set of pdfs 
and a s at each order; the 'NNLO' set relies on an ap- 
proximate set of NNLO splitting kernels • 
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FIG. 1: Fhe CMS rapidity distribution of the lepton pair 
produced in pp collisions at LO (lower band), NLO (mid- 
dle band), and NNLO (upper band), for parameter choices 
relevant for fixed target experiments, along with E866 data. 
The upper (lower) edge of each band denotes the renormal- 
ization/factorization scale choice fi = M/2 (fi — 2M). 

In Fig. 1, we present the center-of-mass system (CMS) 
rapidity distribution of M = 8 GeV lepton pairs pro- 
duced in pp collisions at ^fs = 38.76 GeV, along with 
data from the fixed-target experiment E866/NuSea jg]. 
We set the renormalization and factorization scales both 
equal to \x. The LO, NLO, and NNLO results are shown 
as bands in the figure, indicating the variation of the 
cross sections between the scale choices \i = M/2 and 
fi = 2M. The NLO and NNLO distributions become 
more sharply peaked at central rapidities; this is due pri- 
marily to the evolution of the parton distribution func- 
tions beyond leading order. The significant scale depen- 
dence of the NLO cross section, which reaches nearly 
25% over the interval M/2 < /j, < 2M, is reduced to ap- 
proximately 10% at NNLO. The magnitude of the NNLO 
corrections depends upon the choice of /z; typically, they 
increase the NLO result by approximately 5-15%. We 
note that the NNLO corrections are drastically reduced 
for the scale choice fj, = M/2. The NNLO corrections 
computed in the so-called "soft" approximation, which 
retains only those terms that are singular in the limit 
z —>■ 1, lead to a onnlo which is lower than the result of 
the full calculation by approximately 20% Q. 

The E866 data points are based on measured xf dis- 
tributions |8j, converted from xf to Y with the aid of 
the pj_ distribution [l2j . A ±6.5% normalization error 
common to all data points is not shown. The data lie 
somewhat below the NNLO prediction at smaller Y, and 
rise to meet it at large Y, a trend also visible at other 
values of M. Recall that the antiquark distribution is 
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derived partly from the DY process, fit to the NLO xf 
distribution. Fitting to our NNLO rapidity distribution 
instead may result in a smaller antiquark distribution 
function at moderate x. Further comparison with pd as 
well as pp data is in progress. 
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FIG. 2: The NNLO corrections for the partonic channels qq 
and qg, normalized to the complete NNLO differential cross 
section, for vs = 38.76 GeV, M = 8 GeV, and // = M. 

We now separate our result into its partonic compo- 
nents. The qq and qg pieces contribute the majority of 
the result; the remaining channels are a factor of 50-100 
smaller. The magnitude of the NNLO result is deter- 
mined by a significant cancellation between the qq and 
qg channels. We illustrate this cancellation by plotting 
the NNLO contributions of these channels, together with 
their sum, normalized to the complete NNLO differential 
cross section in Fig. 2. The sum of these channels also 
contributes a much flatter correction to the rapidity dis- 
tribution than either piece individually. The qg channel 
contributes a significant fraction of the complete differ- 
ential cross section; combining both the NLO and NNLO 
qg pieces, we find that they account for about 15% of the 
complete NNLO result for central rapidities, and nearly 
40% for larger (Y > 1) rapidities. This indicates that 
the NNLO rapidity distribution is quite sensitive to the 
gluon content of the colliding protons. 

Finally, we discuss the dependence of the perturba- 
tive K-factors upon rapidity. We define the K-factors as 

follows: ^(N)NLO (F) = (J (N)NLO /(7 LO j and K (2) (y) = 

ct nnlo /ct nlo . We present them in Fig. 3. The signif- 
icant variation of both if NLO (F) and K NNLO (Y) with 
rapidity, a nearly 25% change from Y = to Y = 1, 
illustrates that the LO result provides a poor approxi- 
mation to the shape of the rapidity distribution, as does 
the LO result weighted by the NNLO K-factor computed 



from the inclusive cross section. However, the relative 
flatness of indicates that the NLO result does ac- 
curately predict the shape of the distribution; the NLO 
differential cross section weighted by cr NNLO /o- NLO , the 
ratio of NNLO and NLO inclusive cross sections, is valid 
at these energies to approximately 3-5%. This result ap- 
pears rather promising, since it suggests a simple and 
fairly accurate way of incorporating NNLO corrections 
into NLO Monte Carlo event generators by renormalizing 
with constant if-factors. It remains to be seen, however, 
if the same conclusion is valid more generally. 
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FIG. 3: The K-factors K NLO (Y) = a NLO /a LO , K™ LO (Y) = 
a NNLO /a LO , and K {2 \Y) = a NNLO /a NLO , for (i = M. 

In conclusion, we have described a calculation of the 
NNLO QCD corrections to the rapidity distribution in 
the Drell-Yan process. We have introduced a powerful 
new method for the calculation of differential quanti- 
ties in perturbation theory. Although we have presented 
only a specific example of this technique, it is clear that 
this method is of more general applicability; the relation 
between differential distributions and forward scattering 
amplitudes described above enables the use of multi-loop 
technology for the calculation of a large class of phase 
space integrals. We are confident that this technique will 
be succesfully applied to compute other quantites of phe- 
nomenological interest. 

Acknowledgments: This research is partially 
supported by the DOE under contracts DE-AC03- 
76SF00515 and DE-FG03-94ER-40833 and by the Uni- 
versity of Hawaii startup grant. We thank Mike Leitch 
and Jason Webb for useful discussions. C. A. thanks 
the University of Hawaii for kind hospitality during the 
course of this work. 



[1] S.D. Drell and T.M. Yan, Phys. Rev. Lett. 25 (1970), 
316. 



[2] G. Altarelli el ai, Nucl. Phys. B157, 461 (1979). 



5 



[3] R. Hamberg et al., Nucl. Phys. B 359, 343 (1991) 

[Erratum-ibid. B 644, 403 (2002)]. 
[4] P.J. Rijken and W.L. van Neerven, Phys. Rev.D D51, 44 

(1995). 

[5] E.W.N. Glover, hep-ph/0211412 and references therein. 
[6] C. Anastasiou et al, Nucl. Phys. B646, 220 (2002). 
[7] C. Anastasiou et al, arXiv:hep-ph/0211141 
[8] J.C. Webb et q;., |arXiv:hep^'/ 0302019 



[9] F.V. Tkachov, Phys. Lett. BlOO, 65 (1981); 
K.G. Chetyrkin and F.V. Tkachov, Nucl. Phys. 
B192, 159 (1981). 
[10] T. Gehrmann and E. Remiddi, Nucl. Phys. B 580, 485 
(2000). 

[11] A.D. Marti n et al, Phys. Lett . B531 (2002) 216. 
[12] J.C. Webb, arXiv:hep-ex/0301031 



